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LONG-TERM  GOALS 

We  are  part  of  a  multi-institutional  research  team*  funded  by  the  ONR-sponsored  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program.  The  primary  research  goals  of  the  program  are  to  (1)  examine 
time-dependent  oceanic  radiance  distribution  in  relation  to  dynamic  surface  boundary  layer  (SBL) 
processes;  (2)  construct  a  radiance-based  SBL  model;  (3)  validate  the  model  with  field  observations; 
and  (4)  investigate  the  feasibility  of  inverting  the  model  to  yield  SBL  conditions.  The  goals  of  our 
team  are  to  contribute  innovative  measurements,  analyses  and  models  of  the  sea  surface  roughness  at 
length  scales  as  small  as  a  millimeter.  This  characterization  includes  microscale  and  whitecap  breaking 
waves. 
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OBJECTIVES 

Nonlinear  interfacial  roughness  elements  -  sharp  crested  waves,  breaking  waves  as  well  as  the  foam, 
subsurface  bubbles  and  spray  they  produce,  contribute  substantially  to  the  distortion  of  the  optical 
transmission  through  the  air-sea  interface.  These  common  surface  roughness  features  occur  on  a  wide 
range  of  length  scales,  from  the  dominant  sea  state  down  to  capillary  waves.  Wave  breaking  signatures 
range  from  large  whitecaps  with  their  residual  passive  foam,  down  to  the  ubiquitous  centimeter  scale 
microscale  breakers  that  do  not  entrain  air.  There  is  substantial  complexity  in  the  local  wind-driven  sea 
surface  roughness  microstructure.  Traditional  descriptors  of  sea  surface  roughness  are  scale-integrated 
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statistical  properties,  such  as  significant  wave  height,  mean  squared  slope  (eg.  Cox  and  Munk,  1954) 
and  breaking  probability  (e.g.  Holthuijsen  and  Herbers,  1986).  Subsequently,  spectral  characterisations 
of  wave  height,  slope  and  curvature  have  been  measured,  providing  a  scale  resolution  into  Fourier 
modes  for  these  geometrical  sea  roughness  parameters.  More  recently,  measurements  of  whitecap  crest 
length  spectral  density  (eg.  Phillips  et  al,  2001,  Gemmrich,  2005)  and  microscale  breaker  crest  length 
spectral  density  (eg.  Jessup  and  Phadnis,  2005)  have  been  reported. 

Our  effort  seeks  to  provide  a  more  comprehensive  description  of  the  physical  and  optical  roughness  of 
the  sea  surface.  We  will  achieve  this  by  implementing  a  comprehensive  sea  surface  roughness 
observational  ‘module’  within  the  RADYO  field  program  to  provide  optimal  coverage  of  the 
fundamental  optical  distortion  processes  associated  with  the  air-sea  interface.  Within  our  innovative 
complementary  data  gathering,  analysis  and  modeling  effort,  we  will  pursue  both  spectral  and  phase- 
resolved  perspectives.  These  will  contribute  directly  towards  refining  the  representation  of  surface 
wave  distortion  in  present  air-sea  interfacial  optical  transmission  models. 

APPROACH 

We  will  build  substantially  on  our  accumulated  expertise  in  sea  surface  processes  and  air-sea 
interaction.  We  are  working  within  the  larger  team  (listed  above)  measuring  and  characterizing  the 
surface  roughness.  The  group  plans  to  contribute  the  following  components  to  the  primary  sea  surface 
roughness  data  gathering  effort  in  RaDyO: 

•  polarization  camera  measurements  of  the  sea  surface  slope  topography,  down  to  capillary  wave 
scales,  of  an  approximately  lm  x  lm  patch  of  the  sea  surface  (see  Figure  1),  captured  at  video  rates. 
[Schultz] 

•  co-Iocated  and  synchronous  orthogonal  75  Hz  linear  scanning  laser  altimeter  data  to  provide  spatio- 
temporal  properties  of  the  wave  height  field  (resolved  to  0(0. 5m)  wavelengths)  [Banner,  Morison] 

•  high  resolution  video  imagery  to  record  whitecap  data,  from  two  cameras,  close  range  and  broad 
field  [Gemmrich] 

•  fast  response,  infrared  imagery  to  quantify  properties  of  the  microscale  breakers,  and  surface  layer 
kinematics  and  vorticity  [Zappa] 

•  sonic  anemometer  to  characterize  the  near-surface  wind  speed  and  wind  stress  [Zappa] 

Our  envisaged  data  analysis  effort  will  include:  detailed  analyses  of  the  slope  field  topography;  laser 
altimeter  wave  height  and  large  scale  wave  slope  data;  statistical  distributions  of  whitecap  crest  length 
density  in  different  scale  bands  of  propagation  speed  and  similarly  for  the  microscale  breakers,  as 
functions  of  the  wind  speed/stress  and  the  underlying  dominant  sea  state.  Our  contributions  to  the 
modeling  effort  will  focus  on  using  the  data  to  refine  the  sea  surface  roughness  transfer  function.  This 
comprises  the  representation  of  nonlinearity  and  breaking  surface  wave  effects  including  bubbles, 
passive  foam,  active  whitecap  cover  and  spray,  as  well  as  microscale  breakers. 
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WORK  COMPLETED 


Our  effort  in  FY07  has  been  primarily  in  the  detailed  planning  of  the  suite  of  sea  surface  roughness 
measurements  that  we  will  undertake  during  the  Scripps  Institution  of  Oceanography  (SIO)  Pier 
Experiment  scheduled  for  January  6-28,  2008.  During  FY07  we  refined  our  choices  of  the 
instrumentation  needed  to  make  the  measurements  described  in  the  preceding  section,  and  continued 
work  on  the  analysis  techniques  for  characterizing  the  various  roughness  features.  We  participated  in 
the  FY07  RaDyO  scientific  planning  meetings,  which  were  held  in  Montreal  in  October  2006  and  at 
SIO  in  June  2007. 

We  carried  out  data  processing  and  validation  of  our  scanning  lidar  data  gathered  at  the  USACE  Duck 
pier  field  site  at  Duck,  N.  Carolina.  Two  of  these  instruments,  operating  in  quadrature,  will  be 
deployed  in  the  RaDyO  field  experiments  to  monitor  the  large  scale  wave  geometry  (height  and  slope 
components).  We  also  progressed  with  our  development  of  a  robust  ‘individual  wave’  decomposition 
capability  so  that  local  physical  roughness  elements  can  be  detected  and  characterized  along  with  their 
space-time  phasing,  thereby  overcoming  the  classical  Fourier  spectrum  issue  of  bound  versus  free 
wave  contributions  in  assessing  true  physical  sea  surface  roughness. 

RESULTS 


(i)  Instrumentation  proposed  for  the  SIO  pier  experiment,  January  6-28,  2008. 

Figure  1  below  shows  schematically  the  instrumentation  to  be  deployed  in  this  field  testing  phase. 

2008  SIO  PIER  EXPERIMENT  INSTRUMENTATION  SCHEMATIC 


whitecap 
video  _ 
camera 


polarimetric  camera 


: 


PLAN  VIEW 


orthogonal  scanning  lidars  infrared  camera 


Figure  1.  Schematic  of  instrumentation  packages  to  be  deployed  from  the  northern  swinging  boom 
at  the  end  of  the  Scripps  pier.  The  end  of  the  instrumentation  boom  will  be  about  9m  from  the  edge 
of  the  pier  and  about  10m  above  the  mean  water  level.  The  approximate  field  of  view  of  the  various 
instruments  is  shown.  Another  wide  angle  whitecap  video  camera  is  mounted  well  above  the  pier. 
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Banner/Morison  plan  to  deploy  two  orthogonal  line  scanning  lidars,  synchronized  for  zero  crosstalk. 
These  will  be  positioned  on  the  boom  so  that  their  intersection  point  is  within  the  common  footprint  of 
the  polarimetric  (Schultz),  infrared  (Zappa)  and  visible  (Gemmrich)  imagery  cameras  to  measure 
small-scale  surface  roughness  features  and  breaking  waves. 

Zappa  will  deploy  his  infrared/visible  camera  system  (with  blackbody  target,  a  blackbody  controller,  a 
laser  altimeter).  He  will  also  deploy  his  environmental  monitoring  system  (sonic  anemometer,  a  Licor 
water  vapor  sensor,  a  Vaisala  RH/T/P  probe,  a  motion  package,  a  pyranometer,  and  a  pyrgeometer). 

Gemmrich  will  deploy  2  video  visible  imagery  cameras.  One  camera  will  be  mounted  on  the  main 
boom  next  to  our  other  instrumentation  packages.  The  second  camera  will  be  mounted  higher  up  to 
provide  a  wider  perspective  on  larger  scale  breaking  events. 

Schultz/Corrada-Emmanuel  will  deploy  an  instrument  package  located  on  the  boom  that  includes  a 
polarimetric  camera  imaging  the  very  small-scale  waves,  an  autofocus  controller  for  this  camera,  a 
laser  rangefinder  for  the  autofocus  mechanism,  a  polarimetric  camera  looking  up  at  the  sky  and  a 
motion  package. 

The  individual  data  acquisition  systems  will  be  synchronized  to  GPS  accuracy  so  that  the  various  data 
sets  can  be  interrelated. 

(ii)  Duck  pier  scanning  lidar  breaking  wave  profiles 

Our  scanning  lidar  system  was  field-tested  at  the  Duck  pier  site  over  a  wide  range  of  conditions  where 
the  wind  speed  Uio  ranged  from  light  and  variable,  up  to  16  m/s.  Figure  2  shows  typical  scanning  lidar 
data  over  breaking  waves  under  strong  wind  forcing  conditions. 
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(a) 


(c) 


(e) 


(f> 


Figure  2.  Panels  (a)  through  (f)  show  the  temporal  evolution  of  the  surface  elevation  profile 
measured  by  our  scanning  lidar  of  a  breaking  wave  shoaling  in  7m  water  depth  at  the  Duck  pier 
site.  The  observed  wave  geometry  shows  a  very  steep  leading  face,  with  a  maximum  steepness 
approaching  45°.  Note  the  persistent  feature  developing  on  the  forward  face.  The  x  and  z  axes  show 
distances  in  meters.  The  time  interval  between  the  successive  images  is  0.176  seconds. 


The  lidar  was  deployed  from  a  boom  at  a  height  of  9m  above  the  mean  sea  level,  where  the  water  depth 
was  about  7m.  The  lidar  sensing  spot  on  the  sea  surface  had  a  diameter  of  ~0.2m,  it  made  an  altimetric 
height  determination  every  0.5°  and  the  temporal  scan  rate  was  in  excess  of  50  Hz.  With  its  viewing 
window  of  ±45°  to  the  vertical,  the  maximum  horizontal  data  extent  was  ±  9m  from  nadir. 

We  deployed  this  instrument  for  several  days.  As  anticipated,  the  horizontal  extent  of  the  lidar  return 
increased  as  the  wind  strengthened,  because  of  the  greater  number  of  scatterers  that  provide  the 
specular  return  facets  for  the  lidar  to  make  its  time-of-flight  measurement.  Our  experience  confirmed 
that  this  method  will  provide  useful  data  on  the  height  and  local  directional  slope  of  the  dominant 
waves.  This  information  characterizes  the  background  environment  experienced  by  the  short  wind 
waves  (the  sea  surface  microstructure  roughness).  This  information  will  also  allow  accurate  phasing  of 
the  polarimetric  camera  imagery  of  the  sea  surface  microstructure  with  respect  to  the  underlying 
dominant  wind  waves. 

(iii)  surface  roughness  data  analysis 

As  a  component  of  our  data  analysis  effort,  we  continued  our  effort  to  develop  new  phase-resolved 
descriptions  of  sea  surface  roughness.  In  addition  to  state-of-the-art  Fourier  (e.g.  Elfouhaily  et  ah, 

2003)  and  wavelet  techniques,  we  are  investigating  various  ‘riding  wave’  analyses  (RWA).  This 
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approach  envisages  an  ‘individual  wave’  decomposition  capability  so  that  local  physical  roughness 
elements  can  be  detected  and  characterized  along  with  their  space-time  phasing,  thereby  overcoming 
the  classical  Fourier  spectrum  issue  of  bound  versus  free  wave  contributions  in  assessing  true  physical 
sea  surface  roughness. 

IMPACT/APPLICATIONS 

This  effort  will  provide  a  far  more  detailed  characterization  of  the  wind  driven  air-sea  interface, 
including  wave  breaking  (whitecaps  and  microscale  breaking).  This  is  needed  to  provide  more 
complete  parameterizations  of  these  processes,  which  will  improve  the  accuracy  of  ocean  optical 
radiative  transfer  models  and  trans-interfacial  image  reconstruction  techniques. 

RELATED  PROJECTS 

The  present  project  is  related  generically  to  our  current  ONR  sea  surface  wave  project  in  the  CBLAST 
Hurricane  DRI  entitled:  ‘Wave  breaking  influence  in  a  coupled  model  of  the  atmosphere-ocean  wave 
boundary  layers  under  very  high  wind  conditions’.  While  the  wind  speed  regimes  in  RaDyO  and 
CBLAST  Hurricanes  are  very  different,  common  elements  in  these  two  projects  include  the  need  to 
better  understand  and  parameterize  the  breaking  process  and  how  it  occurs  at  the  different  wave  scales. 

Our  CBLAST  effort  has  resulted  in  a  capability  for  forecasting  wave  breaking  of  the  dominant  waves. 
These  forecasts  validate  well  at  moderate  wind  speeds  of  around  12  m/s.  Validation  of  hurricane 
breaking  waves  awaits  the  data  processing  by  other  CBLAST  Pis.  Our  effort  has  also  highlighted  the 
need  to  better  understand  breaking  at  the  shorter  scales,  where  the  breaker  frequency  statistics  appear 
to  fall  off  towards  shorter  scales.  However,  balancing  wind  input  to  short  waves  with  breaker 
dissipation  rate,  present  modeling  suggests  that  short  wave  breaking  statistics  should  increase  towards 
shorter  scales.  We  will  revisit  this  issue  with  the  new  insights  provided  by  our  RaDyO  datasets. 
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